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SUMMARY 


The small-disturbance equation for flow at Mach number 1 is studied 
with the objective of determining closed analytic solutions representing 
possible fluid motions. Two families of exact solutions appear. The 
first, which represents purely supersonic flow, is, in the two-dimensional 
case, the transonic equivalent of the Frandtl-Meyer expansion and leads 
naturally to simple wave syst ems over two-dimensional surfaces. The com- 
panion solution, in the case of rotational symmetry, yields a conical 
field which coalesces onto the axis. A second family of subsonic solu- 
tions is also calculated. In the case of rotational symmetry, a source- 
like flow with constant mass flux in the downstream half -plane results. 

It is possible to patch these solutions together in such a way as 
to simulate the flow over the rear half of an infinite body. When initial 
conditions on the body dictate sonic speed, expansion waveB first appear 
and accelerate the motion. An abrupt change to the subsonic regime is 
then produced by a single shock of parabolic shape. 


INTRODUCTION 


This paper presents a few exact solutions of the transonic differen- 
tial equation for the perturbation potential when the free-stream Mach 
number is 1 . Such examples, as in the case of Ringleb’s solution (ref. l) 
of the more general potential equation, are of some interest if only by 
virtue of their rarity, and one can hope that they may provide some 
insight into the possibility of introducing further simplification in 
the methods used to solve the basic nonlinear equation. The approach 
used here is effectively a separation of variables technique with the 
additional assumption that the perturbation potential times seme power 
of lateral distance is a function of a single variable that is related 
simply to lateral and longitudinal distance. Since no boundary conditions 
are imposed, it becomes relatively easy to generate solutions and one 
finds that both two-dimensional and axially symmetric flow fields can be 
treated in the same manner . 

The solutions yield, in all cases, either purely supersonic or purely 
subsonic flow. The method of analysis, in fact, thwarts all attempts to 
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attain a mixed type of solution except for the case in which initially 
supersonic flow is joined hy means of a shock wave onto the downstream 
portion of the subsonic solution. The details of the joining are given 
in the final section on applications. 


IMPORTANT SYMBOLS 


a 


C 


P 


k 

M 

P 


u,v,v 


U 

V 

x,y,z 


7 

^max 


P 

T 


4 > 

<P 


a 

h 


speed of sound 

p - p n 

pressure coefficient , 

I 

coefficient^! transonic potential equation (see eq. (lO)) 

Mach number, — 

’ a 

pressure 

perturbation velocity components in x,y,z directions 
total velocity component in free- stream direction 
total velocity 

Cartesian coordinates (if axial symmetry is present, y is radial 
coordinate, see eq. (13)) 

ratio of specific heats, for air 7 = l.h 

maximum angle of flow deflection for attached shock wave 

density 

composite variable (see eq. (15)) 
velocity potential function 
Uq" 1 ® 


Subscripts 

evaluated just ahead of a shock wave 
evaluated just behind a shock wave 
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o evaluated in the free stream 

* evaluated at critical speed (M = 1) 


ANALYSIS 


In this section, analysis "basic to the determination of the exact 
solutions will be presented for two-dimensional and axially symmetric 
flow fields when the free-stream Mach number is 1. The derivation of 
the transonic partial differential equation is first reviewed briefly. 
Particular restrictions on the structure of the fields and the free- 
stream Mach number then lead to consideration of an ordinary differential 
equation that is solvable analytically. The section concludes with a 
discussion of the properties of the derived solutions. 


Basic Equations 


The flow to be studied is, by assumption, produced by slight devia- 
tions in the uniform flow field created by a free stream of velocity U 0 . 
The uniform flow will be directed along the positive x axis of a 
Cartesian coordinate system x,y,z and the perturbation velocity com- 
ponents u,v,w, parallel respectively to the three coordinate axes, are 
assumed small relative to Uo and the speed of Bound a Q in the stream. 
For transonic flow, the free-stream Mach number (M 0 = U 0 /a 0 ) is near 1 
and, in the derivation of the governing potential equation, the non- 
linearity is retained only insofar as the streamwlse component of per- 
turbation velocity affects the result. The continuity equation can thus 
be approximated in the form 

~ [(Uo+u)p(Uo+u,0,0)] + !- (v Po ) + |- (v Po ) = 0 ( 1 ) 

ox ay oz 

where p'(Uo + u,v,w) is density expressed as a function of local velocity 
and p Q = p(U o ,0,0) denotes free-stream conditions. The first term of 
equation (l) represents the gradient of mass flow and because of the 
unidimensionality Is easily evaluated by means of the corresponding 
momentum relation 


= . < 2 > 

where p is p(U o +u,0,0) and p is local pressure. Equation (2), 
together with the isentropic pres sure- density relation p/p Q = (p/p 0 )^ 
yields 
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r , r / \ z i 
* = { 1+ ^ _M ° 2 [ 1 '( 1+ ^)]} 


(3) 


where 7 is the ratio of specific heats. The mass flow is, thus, from 
equation ( 2 ), determined by 


(U 0 + u)p 

U0P0 


( 1 + C*{ 1 +z i iMo l 1 -( 1 <)]} 




W 



and it remains 
wise gradient. 


to evaluate its stream- 


Sketch (a) shows a plot of equa- 
tion (4) when Mo = 0.8. As is to be 
expected, the maximum mass flow occurs 
when the local Mach number is 1. The 
portion of the curve in the neighbor- 
hood of this maximum is the part 
presently of interest. If equation (4) 
is approximated by an expansion 
including only terms of first order 
in u/Uo, the mass flow Is repre- 
, sented by a straight line tangent to 

Sketch (a) the curve at the point u = 0, as 

shown in the sketch. Equation (l) then assumes the form associated with 
linearized compressible-flow theory. A more exact approximation is sought 
here, and to this end the quadratic dependence on u/Uo is retained. 

The mass-flow curve Is thus to be represented by a parabola of the form 


( y u)p - i+U-mo 8 ) f) 

UqP 0 Up C \Po/ 


(5) 


where k is a constant yet to be determined. 

In equation ( 5 ) the ordinate and slope of the parabolic approximation 
correspond to free- stream conditions but some arbitrariness remains in the 
determination of k. The most obvious choice would appear to dictate 
higher-order contact between the two curves, and this leads to the value 

U 0 k = M 0 2 [3- (2-7 )Mo 2 1 (6) 

A more interesting possibility follows if the vertex of the approximating 
parabola Is fixed at M = 1, although the value of the mass flow itself 
may not be exact there. This is the approach used by Oswatitsch (ref. 2) 
and leads to the result 
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U Q k = U 0 


(1-Mp a ) 

a * “ U 0 


( 7 ) 


where a* is the critical speed of sound determined "by M = (U 0 +u*)/a^ = 1. 


A third possibility for the determination of 
known relation 


_ 

Uo " 


1+2 


(1-M 0 £ ) I 

(7 + 1)Mo 2 J 


1/2 


k 


follows from the 

( 8 ) 


If equation ( 7 ) is combined with the first-order expansion of equation (8) 
near Mq = 1, one gets 

Uok=(7+l)Mo 2 (9) 

When equation ( 5 ) is used, together with equation (l) , and a pertur- 
bation velocity potential $(x,y, z) is introduced, the transonic differ- 
ential equation takes the form 

$xx(l ~ ^o 2 - k$x) + ®yy "^®zz = 0 (10) 

where k is given by equation (6), ( 7 ) , or ( 9 ). 


Sketch (b) shows a comparison, at Mo = 1.2, of the three parabolas 



Sketch (b) 


6 


HACA TN 3868 


described above , as well as the more exact curve for mass flow given by the 
isentropic relation of equation (h) . The differences between the curves 
are of the same order of magnitude as those of terms neglected in tran- 
sonic theory and the final choice must be based on additional considera- 
tions. The form of k in equation (9) was proposed originally by Spreiter 
in reference 3 > a 111 ! its superiority in correlating experimental data by 
means of the transonic similarity rules was clearly established at that 
time. Its use will be adopted here for this reason and also because, as 
shown in the following paragraph, it appears as the most natural form to 
relate the transonic shock-polar and mass-flow relations. 

Equation (10) applies to transonic flow in regions for which the 
flow experiences no discontinuities and the gradients of velocity in the 
equation are finite. Transition through a shock wave is governed by a 
finite-difference relation between the components of perturbation velocity 
(u a ,v a ,w a ) ahead of the shock wave and the components (ttjj/Vb/wb) behind 
the shock wave. As given in reference 3 this equation is the transonic 
shock polar 


(l-Mo a )(u a - u-b) 2 ^- (v a - Vb) 2 + (w a - w^) 2 = Mq 2 ( ^p) (u a - Ufa,) 2 

° k^ — - - — ^(ua-u-b) 5 (11) 


where k in the final term is as defined in equation (9) • Consider a 
normal shock wave in the flow with velocity U Q ahead and U 0 + Ub 
behind. Since the mass flow is continuous through the shock wave and 
since in transonic theory the continuous portion of the flow is repre- 
sented always as a potential field with uniform stream conditions, it 
follows that when Ua^Vg^v-b vanish, the value of u^ downstream of 
the shock must correspond to unity on the mass-flow parabola. In 
sketch (b) , unit mass flow is represented by a horizontal line through 
u = 0. The intersection of the line and the parabola corresponding to 
the k of equation (9) is at 


^ = 2 (l " M o 2) 

U o (r+i)Mo 2 

and this agrees with the value of u^ given by equation (ll). This 
result can, in fact, be generalized easily to show that for any normal 
shock the discontinuity predicted by the shock polar agrees with the 
discontinuity given by equation (5) when the mass flow is held fixed. 
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■Exact 




.06t 




04 


EES. 


-20 


A comparison between the tran- 
sonic shock polar of equation (ll) 
and the exact polar is shown in 
sketch (c) for the case Mq = 1.1. 

The inner curves apply to shock 
discontinuities for which free- 
stream conditions hold on the 
upstream face; the two polars are 
thus applicable , in particular, to 
bow waves in a supersonic free 
stream. It is theoretically pos- 
sible to treat shocks at positions 
for which the supersonic free stream 
has been accelerated to a higher Mach 
number associated with the longitudinal perturbation velocity ua. The 
outer curves of the sketch are drawn for the case Mq = 1.1 and ua/Uo = 0.0^. 
It is obvious that Ua/U Q must, in general, be kept very small if the 
shock analysis is to be based on free-stream conditions. 

From equation (ll) the maximum deflection angle accommodated by the 
shock wave can be shown to be 



Sinflx — 


(M a 2 . 


1 ) 


3/2 


3 41 Mo 2 ( 7 +i) 


(12) 


where Mg. is the Mach number on the upstream face of the shock. It 
is clear from an inspection of the curves in sketch (c) that this result 
provides excellent agreement with exact theory when Mq = Ma but becomes 
less accurate as the difference between the Mach numbers increases. 


Particular Solutions When Mq = 1 


Subsequent application will be limited to cases where the flow field 
is either two-dimenBional or has axial symmetry. The differential rela- 
tion for the perturbation potential is, from equation (10), 

(l-M 0 2 -k$ x )$ xx +$ yy +| = 0 (13) 

where a = 0 for two-dimensional flow and cr = 1 for axially symmetric 
flow and y is a Cartesian coordinate in the two-dimensional case and 
the radial coordinate in the case of axial symmetry. When the free- 
stream Mach number is 1, all values of k reduce to 

ko = (7 + D/Uo 
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and by means of the transformation 

$ = U 0 cp 

the differential equation becomes 

- (r+ 1) V 3 cx + V + § V 0 (l4) 

Following Glider ley and Yoshihara (ref. 4) a solution is sought of 
the form 

P-pfW , (k 0 U 0 ) l/3 r = (15) 

If the resulting ordinary differential equation iB to involve only t 
and f ( r) , the condition m + 3 n = 2 must apply. The perturbation velocity 
components then become 


*X xy2-3n 


f, ( T > 


9 y =- 


- ~ - 3 - n 4 f"f(T) + 

y 3(i-n) L 


n 


(2 - 3n) 


Tf*(T) 


(l6a) 

(l6b) 


The ordinary differential equation for f(T) can now be written as 

ftf" = [3(1- n) - a]( 2 - 3n)f + [5n(l- n) - nc]Tf t +n2T2f*' (17a) 

This is a homogeneous differential equation, of degree of homogeneity 3> 
so its order may be reduced by unity (see ref. 5* par. 55) • Aa alternate 
approach is possible, however, since the parameter n is yet to be fixed. 
Rewrite equation (l7a) as 


§ — U’(t)] 2 = n2 JL (T2f)+[5n(l-n) - W 5 - on] A (rf) + 
d ar dT^ at 

[2(l- n)(3 - 7a) + 2n 2 +2c(2n - l)]f (1Tb) 

A first integral of this equation can be written immediately if n is 
chosen so as to make the coefficient of f vanish. Solution of the 
resulting quadratic equation in n yields the roots 
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n 



1=- 1 

k 


For the two-d im ensional case, this results in two possible values for n, 
namely l/2 and 3 /^j "but for the axisymmetric case, where a - 1, the 
roots coalesce at n = l/2. The first integral of equation (l7b) is 


f ,£ = 2n 2 — ( T £ f ) + 2n(5 - 9 n - q) rf + const 


In order to deal with situations in which the flow is undisturbed 
for x < 0, the condition that the streamwise perturbation velocity com- 
ponent qp x vanish at x = 0 will be applied. From, equation (l6a) , this 
implies that f*(0) = 0 also and the value of the constant in the last 
equation becomes zero. The resulting form for the first integral of 
equation ( 17 &) is therefore 


fi2_ 2n 2 T 2ft_ 2n(5-7n-cr)Tf = 0 (l8) 

It will be noted that equation (l8) is again homogeneous, and thus can 
be integrated. It remains then to consider separately the three possible 
cases (cr = 0; n = l/2, 3 /k and a ~ lj n = l/2) and to examine the form 
and properties of the relevant solutions of equation (l8). 


Properties of the Particular Solutions 

Two-dimensional case, n = l/2 .~ Setting a = 0, n = l/2 in equa- 
tion ( 18 ) , one needs to solve the equation 

2f ,a - T2f* - 3Tf = 0 (19) 

and this is readily carried out by means of the substitution given in 
reference 5* First, however, consider the equation written in the form 

f 1 = I ± -7 t 4 + 24rf) ( 20 ) 

It can be seen that f(r) must be an odd function of r, positive when t 
is positive. The ambiguity of sign then leads to two possibilities: if 

the upper sign Is chosen, f 1 > 0, and the choice of the lower sign leads 
to f* < 0. By reference to equation (l 6 a), it follows that these pos- 
sibilities lead respectively to the conclusion that <p > 0 or <p < 0 , 

hence to purely supersonic or purely subsonic perturbation fields, since 
the stream velocity is sonic. These possibilities will be distinguished 
with respect to f by designating them respectively as and f 2 . 
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The solutions of equation (19) are expressible in the formB 


fl(T) , i T . , hill (T 3 + 24f 1 )' L/a + 3 T : 

3 [( T s + 2 lif 1 .) l/s +T s, *] : 


3/2 

1/ 3 


(21a) 


f2( r ) - 2hf 20 


a (t 3 + 2 Uf g ) 1/g + T^/z 
jT3 + 2k-f z ) :L/2 + 3 t 3/2 _ 


(21b) 


where the additional subscript o means that the function in question is 
evaluated at t = 0. A particularly interesting solution of supersonic 
type results when f 10 = 0, namely 

fi(r) =ir3 (22) 


Physical aspects of the flow field determined by this solution will be 
considered later. Supersonic-type solutions given by equation (21a) that 
do not have the property f 10 = 0 will not be considered further here. 

It is apparent, however, that they behave asymptotically like (1/3 ) t 3 , 
the solution to be used. 


Consider, next, the solution (21b) which is of subsonic type, f 2 f < 0. 
The formulation, as given, is implicit but it is not difficult to com- 
pute the course of f 2 (-r) because of the following relation, a consequence 
of the homogeneity. Let 6 = 24f 2 (r) /t 3 ; then equation (21b) can be 
written 



t. ■ ._«? 76 (n-9) 1/a +l 


-20 


9 [(1 + 0) 1/S 


+ 3 


T 


(23) 


so that if values of 9 are assigned, 

T may be computed and hence f 2 , since 


■ £ 


Sketch (d) shows the behavior of 
f £ (r) for the initial value f ao * 1« 
Because of the homogeneity of f (t) 
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in r, it is only necessary to confute * f 2 (r) ^ or orLe initial value T20J 
other cases then follow from the transformations t=ct*, f g ( T ) = c 3 f 2 *(r*). 

For later analysis, it is convenient to have simple expressions for 
the "behavior of f 2 ( r) for smal 1 and for large values of ■ the argument r. 
Such expressions are readily calculated by use of the defining equa- 
tion (19). The required formula for small and large t are, respectively. 



^2(1-) 


1 

k 


-20 


1 27 r 4 

•75-16 f ZO 



(2 to) 


Two-dimensional case, n = 3 /^ «~ The function f(r) is now replaced 

by g(r) in order to distinguish the present functions from those for 
n = l/ 2 . The defining equation is 


8g' 2 - 9^* + 3 Tg = 0 
Solving for g* , it follows that 

S* = ^ [ 9 t*± ( 81 t 4 - 96 rg) 2jz ] 

Thus g( t) must be negative, or such that 

,81 3 
g C —T T 3 
96 


(25a) 


( 25 b) 


The solutions associated with the two signs in equation (25b) will again 
be distinguished by the subscripts 1 and 2, corresponding to the positive 
and negative signs, respectively. These solutions are 


8 i(t) = -| t 3 


(96g 


10 


1 IT (81t 3 - 96gir ) 1/2 + 7t 3/2 

[(81t 3 - 96Tg i ) l/2 + 9 t 3/2 J 


( 26 a) 


g2( T ) 


/ S go 8 V /9 ( 81 t 3 - 96Tg 2 ) 1/2 + 9T 3/g . 

[(8lT3-96Tg E ) 1/2 - + 7T 3/2 ] 7,S 


(26b) 
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-.8 


- 1.2 


- 1.6 


- 2.0 


-2.4 l 


gi(r) 


/82(f) 


Sketch (e) 

feM » 8 *,- (^f 2 ) 


1 / 2 


It is of interest to remark 
that the supersonic-type solution 
for which g 1Q = 0 coincides with 
that found in the previous case. 
The solution of subsonic type is 
shown in sketch (e)j it does not 
vanish for large t as did f2(r) 
hut, in fact, grows as t 1 ^ 3 . 
Approximate formulas for small and 
large t are again easy to find, 
being 


T 3 / 2 + -X- T 3 + 

48 


(27a) 


/ d\l/9 A, o\ 2/9 

- ( T ) ~_£ (§SS?) W 3 --l ($22?) -k 

8zK) 4V192 J 16 V 1 92/ 


+ • 


(27b) 


Axisymmetric case.- For this case the functional symbol F(t) is 
used to replace f (t) ." The defining differential equation is 


2 F * 2 - t?F * - tF = 0 

which can be rewritten as 

F* = X [t 2 ± (t 4 + 8tF) i/2 ] 

Solutions of equation (29) are 


Fi(t) = 9 t3 + 


1 Aio* 
3 V 8 , 


1/5 3 (T 3 + 8F 1 ) l/g +9T 3/2 

[( tS + SFx ) 172 -^^ 2 ] 375 


(28) 


(29) 


(30a) 
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FeCx) - 6 M 1 ' 3 

[3(t3 + 8f 2 ) i/2 + 5 t 3/2 ] 


(30b) 


Once again a simple, supersonic- 12 
type solution results when F 10 = 0 , 


-! t3 


(31) 


The subsonic solution can be calcu- - 
lated by the same device used pre- 
viously and has the form shown in 
sketch (f) . Expansions for small as 0 ( 
well as large values of t are 


s \l/2 f \1/ 2 

Fit(x) -Fso-f (-f 2 ) t3/2 + ^ ’ S/a+ - • • (32a) 



Sketch (f) 
l/ 2 


F 2 (t) (SfFso 4 ) 173 (SFso 8 } 173 ^ 


(32b) 


The formal similarity between the two- and three-dimensional cases 
is particularly striking for the simple, supersonic -type solutions since 
they differ only in the magnitude of the coefficient of t 3 . It will be 
seen in the next section, however, that the two-dimensional solution is 
of more general interest since It Is directly related to the study of 
expansion waves in transonic flow. 


SPECIAL APPLICATIONS 


In this section streamlines and pressures associated with the derived 
solutions will be studied in some detail. In all cases the flow in the 
negative half -plane (x < 0) will be assumed uniform with parallel stream- 
lines. The disturbance fields then produce deviations in slope of the 
strea m l in es In the positive half -plane (x > 0 ) and the flow can thus be 
Identified with modifications in shape on the rear of a semi-infinite 
cylinder or airfoil. 
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Two-Dimensional Flow 


Super sonic-type solution .- Consider first the supersonic solution 
determined "by the function fi(r) given in equation (22): fi(T) = t 3 /3. 

As mentioned previously, this solution holds when = for either 

n = 1/2 or n = 3/4. According to equations ( 15 ) and (l6) , the perturba- 
tion potential and velocity components are 


cp(x,y) 


1 x 3 ! 1 

3y 1/2 (Uoiso)y 3/a ~ 3(7 + 1) y 2 


(33a) 


9 x (x,y) 


1 

(7 + 1) 





q>y(x,y) 


-2 

3(7 + 1) 



(33*) 


Consistent with the approximations made in establishing equation (10), 
the slope of the streamlines is given by 


§2. = 
dx 


^y 


(34) 


and, after integration, this yields 


y* + 


3(7 + 1) 


x 4 = y c 4 


(35) 


where y Q is the starting value of y. The flow therefore expands as 
it becomes supersonic. The solution given in equations (33) is clearly 
of the sort called "simple wave" since it lies adjacent to an undisturbed 
region. The characteristics of the field should then consist of at least 
one set of straight lines, and this is seen to be the case; for, from 
equation (l4) with c = 0, the slopes of the families of characteristics 
are 


= 11 

** [ (7 + i)cp x ] 


and, using equation (33*) > 


& = ± £ 
dx x 


(36) 

(37) 


The characteristics are then 
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y/x = const , xy = const ( 38 ) 

and the first of these families is the required set of straight ' lines . 

The "best known simple wave solution is that of the Prandtl-Meyer 
expansion and, in fact, the present solution is the transonic approxi- 
mation to this exact solution, applying for small angles of flow deflec- 
tion. The accompanying sketch (g) shows a comparison of the two flows 
as to streamlines and second-family characteristics (the first-family 
characteristics are identical, "being radial lines from the origin.) 



Sketch (g) 


The solid l in es are the exact results and the dashed lines represent the 
present approximation . Roughly speaking, agreement is good up to a fan 
angle of about 45°, corresponding to a deflection, or expansion, through 
about 6°. Another comparison that can be made is that of the pressure 
coefficients in the two solutions. The exact result is (for Mq = l) 


Cp 




( 39 ) 


where 7 is the ratio of specific heats, A 2 = (7 - l) / (7 + 1) and 
tan w = x/y so that w is the polar angle measured from the y axis. 
For mrifi.1 1 angles of deflection a, it is easily found that the relation 
between a and w is 
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a 


3(7 + 1) 


CO ' 1 


If the right side of equation (39) is expanded for small values of 10 
and the lowest-order term retained, one gets 


C 


P 


« -=£- u2 
7 + 1 


In terms of a, the last relation yields for pressure coefficient 


„ „ -o, 2/3 (l8) l/3 

° P ~ (7 + 1) 1/3 


m 


In the present transonic small-perturbation theory, pressure coefficient 
is — 


Cp ‘ ' ap * = (r’+x) (*) ‘ Cr + x) “* 

The deflection, <x, of the flow is 

2 /x\ 3 2 

a = = 3(7 + 1) V/ = 3(7 +1) “ 3 


and these last two relations lead directly to the approximation for 
Prandtl-Meyer flow given in equation (4o). This expression is, moreover", 
in accordance at Mq = 1 with the transonic similarity rule for pressure 



coefficient as given, for example, 
in reference 3. 

The nature of the Prandtl-Meyer 
relation between pressure and local 
flow. deflection leads naturally to 
application to two-dimensional expan- 
sion over a curved airfoil surface 
for arbitrary free-Btream Mach numbers. 
The problem resolves itself into the 
determination, from a knowledge of the 
expansion fan for Mo = 1, of a rela- 
tion for pressure coefficient, analo- 
gous to equation (40), for the case 
In which a flow at a Mach number M 0 
Is turned through a small angle. "To 
determine such a formula, one starts 
with the previous solution and follows 
it until the Mach number in the expan- 
sion becomes Mo, at some angle Wo, 
as shown in sketch (h) . The velocity 


Sketch (h) 
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is constant all along the ray at angle 0 ) o to the y axis , and the 
direction of the flow is parallel to the ray at angle Oq to the x 
axis. Introduce now a new potential function i|r(x,y) defined by 


*(x,y) 


(Mo a -l)x + qp 
(r+l)Mo 2 Mo 2 


m 


This function satisfies the transonic equation ( 13 ) for two-dimensional 
flow at free-stream Mach number Mo, and so applies to the flow in the 
present problem as it exists after turning through the angle a Q * The 
boundary conditions are 


= OlQ 2 - (Mp 2 - 1) 
(y + l)Mo 2 


at o)q 


and since w 0 2 = Mq 2 - 1 in the present approximation, one has 


t x = 0 at to 0 


In general, the streamwise component of perturbation velocity is 


w 2 - w Q 2 
(7+1)Mo 2 


The deflection angle is 


-2 


«*>*■ 


* y 3(7+1) M 0 ‘ 


— -(cto + cc.) 


and the initial boundary condition yields 


a-o = t T 

J C0=t)J o 


-2lo q 3 = -2(Mq 2 - 1 ) 3/2 
3(7 + l)M 0 2 3(7 + DMo 2 


Thus the incremental pressure coefficient attributable to the deflection 
from Oq to a Q +a is 


C P 




(Mp 2 -1) 

(7 + 1)Mq 2 \ u o 2 


or, if the relation between tu and a is used. 
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Cp = 


(Mo g -i) f r 3(7+ i)Mq 2 
Cr + l)Mo 2 l L 2(Mo a -l) 3/2 


- 2 / 3 -\ 

aj ) 


m 


Equation (42) is consistent -with, the transonic similarity rule for 
pressure coefficient -when k = (7+ l)Mo 2 /Uo is the haBic transonic param- 
eter. It is of interest to remark that the same problem has been treated 
by Lighthill in his discussion of simple wave theory of airfoil flow 
(ref. 6, p. 387)* His result , which was derived as a summation of series 
representations for pressure coefficient in terms of flew deflection angle, 
was 


C 


P 


- 2 


(Mp g - 1) 

(7 + I) 



3(7 + 1) 


3/2 


2 (Mo 2 - 1) 



(43) 


Here, a is of opposite sign than that used in equation (42). It Is clear 
that equation (43) is in a form consistent with the transonic similarity 
rules based upon the value (7+l)/U 0 of the transonic parameter k. 


Sub sonic -type solutions .- Consider the two-dimensional solution for 
the case n = l/2. The differential equation defining the streamlines 
of the field is 

£"'Py“t3^t f a( T > + rt a’( T >> (**■») 


This equation must be solved numerically and for the reason of avoiding 
continual conversion between x, y, and t. It .is convenient to change the 
independent variable from x to t, This gives 


= _2y f at Tgg 1 

dT 4(U 0 k 0 )" 1 ' V + r(f 2 + rf 2 » ) 


(44b) 



A plot of the quantity (f 2 + Tfo’ ) 
versus r is shown in sketch (i) . 
The quantity vanishes at one value 
of r, say t 1# and this is easily 
found to be at the point where 
f2 ( T i) = T i 3 (hy using eq. (20) for 
f 2 * in terms of f 2 ) or, using 
equation (21b) , where 



^20 


2 


Sketch (i) 
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Tims, the streamlines have a stationary point at the zero of f 2 +Tf’ 
and, in the x,y plane, a constant value of t becomes a parabola 

x 2 = [(Uoko ) 17 3 tJ y 

Sketch Q) shows some results of integrating equation (44b), and 
retransforming the Independent variable from t to x, as well as the parab- 
ola on which the minimum points of the streamlines lie. In addition one 
can derive various other relations for the streamlines as approximations . 



Sketch (j) 

For example, for very small values of t (or x) , one can use the first 
term only of the series (24a) for f 2 ( t) , and derive the starting form 
of the curves. This gives 

y 572 = y 0 5/2 “ \ f so* 

where 


y Q = y(°) 


On the other hand, to determine the asymptotic form of the streamlines, 
the first term of the series (24b) can be used together with equation (44a). 
Thus 


y 0 


.. _ 3(Uoko) 2 1 

,-y — § — 20 
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where y TO is the ultimate height of the streamline . In this particular 
example , the streamlines approach their asymptotes from "below. Finally, 
for large starting values of the streamlines , y Q » 1, a control surface 
method would "be applicable for the calculation, since the variation dies 
out rapidly with increasing y o; and the streamlines are nearly straight. 

It remains to calculate the pressure distribution along a streamline 
in the flow. The expression for pressure coefficient is given explicitly 
as r 

gp-.apx- ~?1 n CW 

(tfc£o) y 



Sketch (k) 


The results of the numerical solu- 
tion of the differential equa- 
tion (44b) for the Btreamlines are 
useful for this pressure calculation, 
since therein y is known as a 
function of t. Results for the 
streamlines having y p = 1, based 
on the function f 2 (Tj for which 
f ao =1, are shown in sketch (k) . 

It is not difficult to determine 
that near x = 0, 

EZ-l. 

dx \/x 


so that the pressure gradient is infinite at the y axis, on every 
streamline. This effect Is attributable to the finite deflection angle 
on each streamline at x = 0. 

The other two-dimensional solution: (n = 3/4) can be discussed in 
much the same way as for the preceding case (n = l/2) . The differential 
equation for the streamlines is 


Again, the independent 
in 

: 

dT 


& = ■ 
dx 3*ry 


3 / 4 


S2 ,2 (t) 


variable can be changed from x to t, resulting 


4 2 ga >2 ( T ) 

3 T 2(U 0 ko) l/3 y + g 2 ,2 (r) 


f 
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Re stilt b for numerical solution of 
this equation, for several initial 
values of the lateral coordinate 
y Q , are indicated in sketch (z). 

The initial value of g 2 (-r) was 
S20 = -I* It is noticed that there 
are no stationary points on these 
streamlines, but the corner at 
x = 0 persists. Pressure on the 
streamline for which y Q = 1 is 
shown in sketch (m) , and the singu- 
larity in the pressure gradient at 
x = 0 still exists. 

Combination of subsonic and 
supersonic solutions .- It is of 
interest to consider the possibility 
of combining the previous solutions 
through the mechanism of a shock 
wave. Thus, consider again the 
undisturbed sonic flow in the half- 
plane x < 0, and let it expand to 
become supersonic for x > 0 . If 
this supersonic flow Is followed for 
some small distance, it can be shown - 
that a shock wave of determined 
strength can be specified in such a .8 l 
way that the downstream field corre- 
sponds to known results. Assume now 
that the simple-wave, supersonic solution is to be followed by the sub- 
sonic solution for the case n = l/ 2 . The discontinuities in the flow 
must satisfy the two-dimensional fora of the shock-wave relation given 
in equation (ll). Since Mq = 1 , the condition can be written in the 
form 



Sketch (z) 



Sketch (m) 


(’Va-’VbJ ‘ (V%) 


( 46 ) 


In terms of the auxiliary functions fi(T) and f 2 (r), equation (k 6 ) 
becomes 

(f 1 - f 2 + Tf x » - rf 2 * ) 2 - 2(f !* + f 2 l ) (f JL* - f 2 ’ ) Z - 0 
which factors to give 


(f x - f 2 ) (f ! - Tf 3.' - f 2 +Tf 2 «) = 0 
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The shock-wave relation is thus satisfied if 

fi(T) = f 2 (r) . (47) 

(The other factor cannot vanish since f ]_ - rf j/ is a negative quantity 
and f 2 - Tf 2 * is always positive.) Since equation (47) will hold for 
some fixed value of r, say r 8 , it is clear that the shock wave in the 
physical plane must have a parabolic shape, with the equation 

x 2 = (Uok 0 ) a/ 3 T s a y (48) 

Further, it is not difficult to determine, by setting f 2 ( T ) = t 3 /3 in 
equation ( 21 b), that 



In order to see if the condition represented by equation (47) is indeed 
sufficient for the existence of a shock wave, velocities can be calculated 
on both sides of the parabola of equation (48) in directions tangent to 
and normal to the .parabola. One then finds that the tangential component b 
are continuous across the parabola and that the normal components fulfill, 
to the order of transonic small-perturbation theory, the condition 


where V n and represent components of the total fluid velocity normal 
and tangential to the shock wave, respectively, and the subscripts a 



and b refer, to conditions immediately 
upstream and downstream of the shock 
wave. 

A final condition to be met, if 
the joining of the two flow fields 
by a shock wave is to be meaningful, 
is that the shock wave be attached 
to the wall. Equation ( 12 ) gives an 
expression for the shock-attachment 
angle &max> sketch (n) shows a 
comparison of the exact values of 
Smny as derived from oblique -shock- 
wave theory with the approximate 
values given by equation ( 12 ) with 
Mo “ 1 as a function of the local 
Mach number just ahead of the shock 
wave. Mg,. This sketch re-emphasizes 
the necessity for confining the 


2 * 
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application of this perturbation analysis, where shock conditions are 
referred to free-stream quantities, to very small values of the perturbed 
quantities. Now, to insure that deflection angles obtained in matching 
the supersonic and subsonic flow fields are admissible, it is necessary 
to determine the angle at which streamlines o£ both parts intersect at 
the shock wave. If A is the angle of inclination of a streamline, its 
slope is 


tan A = & = Vy = 

dx y - 2 y 3/ 2 


The discontinuity across the shock wave at any point thereof is then, 
making use of equation ( 20 ) and the shock condition (47) , 


ian(Ab ~ ^a) ,3/2 T e 

4y 1 __l£_ 

l 8 y 3 

. (3 1/a /a)(Wy 3/a > (50) 

1 - (2f 20 2 /27y 3 ) 

Notice that the deflection decreases with the 3/2 power of distance from 
the x axis. The angle determined by equation ( 50 ) is to be compared 
with 5m»Y, using as the local Mach number in equation (12), or in 

sketch (n) , that attained at the end of the supersonic portion of the 
desired combined stream l ine. The local Mach number M is determined 
from the transonic relation 


1 - M 2 « 1 - Mq 2 - (Uok)cp x (51) 

which appro ximat es the coefficient (l-M 2 ) of cp^ in the exact potential 
equation by the coefficient of cp^ in equation ( 13 ) (where U 0 cp = $). 

The calculation of an illustrative example is most easily carried 
out by stenting first with one of the subsonic streamlines already cal- 
culated and matching to it, at the point of intersection with the shock 
wave, a supersonic streamline. The latter is readily determined from 
the differential equation (34) when one uses as the boundary condition the 
values of x,y obtained from the intersection of the wave and the 
subsonic streamline. 

In order not to extend the supersonic flow to angles of deflection 
that exceed the range of validity of the theory, it Is necessary to 
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choose f ao properly. The results of calculations for which f ao = l/8 
are shown in sketch (o) . The subsonic streamline was taken as the one 



x 


Sketch (o) 

for which y Q = 1. The shock wave is a portion of the parabola 

x 2 = (U 0 k 0 ) 2|/ 3 T s 2 y 

which becomes, after evaluating r s from equation (49) , 

x 2 = 0.493y 

From equation (50) , the calculation to determine change in deflection 
angle at the shock wave is as follows 

tan(A b - A a ) = 0.111 , A b - A a = 6.3° 

In order to calculate the local Mach number.. Ma at the end of the expan- 
sion region of this flow, we use the relation (51) which, for Mq = 1, 
gives 

Ma - [l+Cr+UtPsa ] 1 ' 2 

and, from equation (33h) , 

1/2 


/ x B \ z ~ 

* L l+ (?5) J 
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With the value of Mg, thus determined, sketch (n) can he entered, and 
it is found that the angle - A a just determined lies above the shock- 
attachment limit, so that the shock wave is attached to the wall. 

In connection with the wall shape and shock wave of sketch (o) , it 
should be mentioned that the analysis leads to the so-called "strong 
shock" and that it is not possible, in this two-dimensional case, to 
derive a shape characterized by a "weak shock." This result can be found 
by a consideration of the perturbation shock relation (ll) . The region 
of the shock polar that lies between the value of cp x ^ corresponding -to 

maximum deflection angle and cp v =0 corresponds to the weak shock 
region* It is easy to determine that in this region 

- — < % < 0 

it can be shown that this condition corresponds, 
t, to 

■*=>(§^0) >V 

where r s is given in equation ( 49 ). The parabola in the xy plane for 

which cp-^ > - (Txa/3) therefore lies out Bide the shock-wave parabola 

determined by t b . Thus, the weak-shock condition can never be achieved 
in this case. 

On the supersonic portion of the streamline just illustrated, pressure 
coefficient is, from equation (33b) 


Now, using equation (16) , 
in terms of the variable 


°P - -2cp x 


2 

7 + 1 



where y is the height of the streamline, y = y(x) . If the point where 
the shock wave intersects the streamline is denoted (x s ,y s ), pressure 
coefficient on the upstream side is 


Cp -■ 


7 + 1 


y s 4 + 


3(7 + 1) 


&S 4 - X 4 ) 


-ii/ a 


( 52 ) 
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In the Bubsonic portion of the field, 
calculation of pressure is not so 
simple, but is given by 



2.0 


C-n - 


(7 + I) 


U 3 .. 


f 2 ’( , 0 


Sketch (p) 


where y = y(x) is the calculated 
height of the subsonic streamline. - 
Results of the pressure calculations 
on the wall shown in sketch (o) are 
given in sketch (p) . 

It is not possible to combine so 
directly the subsonic and supersonic 
flows when the other two-dimensional 
solution, g( t) , iB to be used. The 
reason is that since the condition 
for a shock wave is again g x = g 2 , 
the simple super sonic -type solution 
g x = t 3 /3 is not applicable since 
the subsonic solution g 2 is always 
negative and no intersection occurs. 
Further application of these results 
would involve detailed examination of 
the more general supersonic solution 
given by equation ( 26a) . 


Axially Symmetric Flow 


Supersonic solution .- A single solution exists in the axially sym- 
metric case, corresponding to the value l /2 for the index n. There are, 
of course, still the solutions of supersonic type and of subsonic type, 
as given in equations (30) and (31). A simple supersonic solution exists 
again, namely, from equation (31) 


Fx(t) - I t3 

The corresponding potential and velocity functions are 


9 


■ -ki F ^ T > 


<p = 


2 x 3 
9(U 0 ko) 7 Z 


9 * = 


x 3(Uoko) V> 


- 9(Uoko) 


( 53 ) 
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so that certain aspects of this flow" are similar to those in the two- 
dimensional supersonic case previously considered. A definite difference 
does exist, however , in the fact that the two-dimensional stre am l in es 
could he carried to the limit as y approaches zero and resulting in 
this case in a Prandtl-Meyer corner flow. This limiting process fails 
three-dimensionally since the flow turns inward on itself. Finite values 
of streamline radius are therefore required and the flow must be modified 
downstream by a shock wave and a readjustment of the streamlines. 


Subsonic solution .- The subsonic solution for the axially symmetric 
transonic problem depends on the function F 2 ( t) defined in equation (30b) 
and shown in sketch (f). Streamlines corresponding to this case can 
again be calculated by solving the differential equation 

1 [ F ( T ) + tF 2 » ( t) ] 

dx 23^ 2 ^ 

By means of the differential equation ( 28 ), this can be written 

& — ( - ^j - 3 5V a (T) (5te) 


Again, it is convenient to change 
independent variables from x to r, 
and the resulting equation is 


djr = _ 2 y F 2 ' 2 (t) 

dT T 2 (U 0 ka) _ 1 / 3 y 2 + F 2 , 2 (t) 

(5 ^b) 


In this case, there are no station- 
ary points on the streamlines. A 
few results of calculations are 
shown in sketch (q) (where ^20 
has been taken equal to l) . These 
streamlines approach their asymp- 
totic height from above, as in the 
two-dimensional case with n = 3 1 ^> 
rather than from below, as in the 
two-dimensional case with n = l/ 2 . 
Pressure calculations proceed as 
previously, the formula being 


°p = “ 


( 7 + 1 ) 


1/3 


F 2 1 (t) 



Sketch (q) 



In sketch (r) , the pressure on the 
streamline for which y Q = 1 is 
shown. 


Sketch (r) 
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Combined supersonic and subsonic solutions .- For the axially sym- 
metric problem, it is again possible to use the simplest supersonic 
function Fi(t; = (2/9)t 3 in con junction. 'with the subsonic solution 
F 2 (t), as given by equation (30b), to produce a flow in which a shock 
wave is imbedded. The condition for existence of a shock wave in the 
axially symmetric case is, from equation (ll) 

(% - %) - (% - %) 

where y now denotes the radial coordinate.. In terms of the auxiliary 
function F(t), this last condition can be written 

(f 1 -f 2 )(f 1 ’ 2 - f 2 > 2 ) = 0 (55) 


Again it can be shown that the first of these factors leads to a proper 
relation of tangential and normal velocities across the shock wave, while 
the second does not. The condition is, therefore, that the functions 
F }_( t) and F 2 (t) be equal at the shock wave. The value of t at which 
this occurs is 

-f(^f <*■> 


and the shock wave is then given in the physical plane by 

x2 = [(Uoko) 373 "^] y 


(56b) 


To illustrate the sort of body that one gets in this case, consider a 
streamline (subsonic) for which y G = 1.2, and computed on the basis 
F 2o = l/8. The procedure. is the same as it was previously, namely to 
determine the intersection, say (xg,y B ) , of this streamline with the 
parabola (eq. (56b)) representing the shock wave, and to use this point 
as a boundary condition for determining the proper solution to the 
supersonic streamline differential equation 

& = — it — (£) 

dx 9(Uoko) W 


(57) 
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Sketch (s) 


The result of such a calculation, -with the parameters y D and. F 20 taken 
as mentioned, gives the figure shown in sketch (s). The angle between 
streamlines on the two sides of the shock wave is, in this case 

*b- « 5-7° 


In this axially symmetric case one gets by using equation (51) with 

Mo = 1, 



and calculation of Ma. shows that the shock wave is attached ih this 
case. In contrast to the two-d im ensional case treated above, it is found 
that in the present axially symmetric case it iB necessary that the shock 
wave be of the weak family, because the critical parabola is determined by 

74 ■ H ’Mf Feo ) < t ° 4 

where t s is given in equation (5 6), 
and therefore lies inside the shock 
wave. Since in both the two- 
dimensional case and the present 
one, the line of demarcation between Cp 
strong and weak Bhocks is distinct 
from the shock wave itself, one 
cannot determine configurations 
having shock waves that are right 
at the point of detachment. 1 

The pressure coefficient on the 
streamlines of sketch (s) is shown 
in sketch (t) . 



Sketch (t) 
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Mass-flow property of subsonic solution with axial symmetry .- All 
of the solutions given here were derived hy means of formal manipulation 
of the transonic differential equation but in the case of the supersonic 
solutions f i(t) = t 3/3, g^T) = t3/3, and F^r) = (2/9)1^ it was easy 
to identify these results with simple expansions from uniform stream, 
conditions. It is evident that the latter solutions could have been 
derived from any one of several methods of attack since their existence 
and nature are intuitively so obvious. The subsonic solutions , an. the 
other hand, are somewhat more complex and. appear to arise rather arti- 
ficially as a result of abrupt deflections imposed on the streamlines 
across the entire plane of flow. In the case of axial symmetry, however, 
the subsonic solution has properties reminiscent of the fundamental solu- 
tion identified with the source-sink potential of linearized supersonic 
theory. From equations (15) and (l 6) the perturbation potential and 
velocity components can be rewritten as 

9 = -ijE F a( T ) (58) 

y 


u - 9 , 


*V(t ) 

y(r + i) 


1/3 


v = q> y 


1 _d_ 

2yS/ 2 &T 


[tF 2 (t)] 


(59) 


From equations (32), limiting conditions yield: 
For t « 1 


u « 


(7+1) 


1/3 


F 20 

2 J 


ft 


1/2 


ri/a 

y 


. 1 (issY* til + 

. (7 + 1) 1 ' 3 W y 374 


V 


1 

2y3/2 


+ . 


For t » 1 


u 



(SEW 1 ) 1 ' 3 


1 

r2y 


+ 


= (9F 


) 1/3 ^ + 


20 


1_ 

x 2 


6 k C3Fao 


1/3 


T 5y' 


3/2 


-g (jho 8 ) 


1/ 3 


X 5 


It is apparent that the origin is an essential singularity of the 
flow field, the rate of growth of the velocity components being determined 
in part by the particular parabola t « const along which one approaches 
the origin. Away from the origin and On the lateral axis u is zero and 
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v is finite but attenuating (as l/y 3 ^ 2 ) with distance. On the longi- 
tudinal axis v is zero and u is finite but attenuating (as l/x 2 ) 
with distance. The flow thus appears to flow into the origin and to 
distort the field in the portion of the plane downstream of the Mach 
line at x = 0 . 


It remains to calculate the rate of mass flow that is associated 
with the solution; if this rate is constant, the analogy with the source 
is more closely established. To this end, consider a cylinder of fixed 
radius Y and extending from the plane x = 0 to x = ». The rate of 
flow is then given by the integral 


R = 


lim 
Z — ^ 00 


Po 


J SnTtvly^dx 

o 


lim 
Z -S* 



a[TF 2 (r)] /_x_\ 

At 


Since d(x/Y 1/2 ) = (7 + 1) 3 dr, the integration is direct and yields 
R = - 2 * Po ( 7 + 1 ) 1/3 lim ZF 2 (Z) =-| rtp 0 ( 7 + l) ^ 3 (9F 2D 4 ) X/ 3 

Z CO 4 

The rate of mass flow is therefore a negative constant independent 
of Y and the solution behaves macroscopically like a sink. As seen in 
sketch (q), the streamlines of the field are drawn toward the longitudinal 
axis, consistent with this concept. 


Ames Aeronautical Laboratory 

Rational Advisory Committee for Aeronautics 
Moffett Field, Calif., Aug. 1 6 , 1956 
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